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SUUM&HT 

Transient  response  data  obtained  from  the  flight  tests  of 
an  F-oQA  airplane  v;ere  analysed  for  the  longitudinal  stability 
derivatives  by  the  Equations  o£  Hot ion  Efethod  in  order  to  de- 
termine whether  or  not  the  method  is  practicable.  Other  con- 
siderations include  the  possibility  of  determining  the  moment 
of  inertia  of  the  aircraft  from  transient  response  data  and  a 
relative  evaluation  of  the  type  of  input  that  should  be  used  to 
obtain  the  transient  data. 

All  data  considered  vrere  obtained  at  a  constant  center  of 
gravity  position  and  essentially  the  same  altitude  and  Ilach 
Tromber. 

The  method  is  considered  to  be  a  practical  means  of  obtain- 
ing the  longitudinal  stability  derivatives  of  an  airplane  pro- 
vided  certain  precautions  arc  observed.  The  accuracy  of  the 
moment  of  inertia  determination  is  somewhat  less  than  that  ob- 
tainable from  ground  tests.  In  general  the  step  function  data 
yielded  better  results  than  the  impulse  data. 

The  values  of  the  stability  derivatives  and  the  inertia 
parameter,  h,  are  as  follows: 


Otep  Input     .impulse  Input 
CL  ,  per  radian       6.22  £.91 


0mr  >  Per  radian      -:,.?6  -0.36 

0 


Gv^     ,  per  radian  — .  .32U  -0.35>7 


C^     ,  per  radian  -0.023  -O.033 

Cmdtf  >  P«P  radian  -0.01(1;  -0.027 

h  +  0.015  -f  0.010 


IlITnODUCTIOU 

In  recont  years  a  groat  deal  of  attention  has  been  directed 
toward  determining  the  stability  derivatives  of  missiles  and 
airplanes.  Steady-state  testing  techniques  and  trind  tunnel 
information  are  used  effectively  to  determine  the  stability 
characteristics  for  subsonic  aircraft,  Hovirever,  with  the  coming 
of  jet  and  rocket  powered  aircraft  and  with  the  resulting  transonic 
and  supersonic  speeds,  it  has  become  increasingly  difficult  to 
obtain  reliable  information  from  v/ind  tunnel  tests*  Also  the 
inherent  difficulty  in  achieving  steady  flight  conditions  in 
the  transonic  and  supersonic  speed  ranges, combined  with  the 
relatively  brief  testing  period  available,  severely  limits  the 
use  of  steady-state  testing  procedures.  It  is,  therefore,  highly 
desirable  that  sone  satisfactory  means  of  analysis  of  transient 
response  data  be  developed. 

Present  day  methods  of  analysing  transient  data  can  be 
broadly  classified  as  follows: 

(1)  Tethods  dependent  upon  the  solution  of  the  equa- 
tions of  motion  of  the  airplane 

(2)  Curve  fitting  methods 

Outlines  of  the  variations  in  each  of  the  above  classifica- 
tions may  be  found  in  Rof.  1  and  2. 

The  primary  objective  of  the  present  investigation  is  to 
determine  the  feasibility  of  analyzing  transient  res "onse  data 
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by  the  "Equations  of  Lfotion  Ifethod"  in  order  to  determine  the 
longitudinal  stability  derivatives  of  an  airplane.  Secondary 
considerations  include  (1)  the  relative  merit  of  the  response 
to  an  impulse  input  as  compared  to  that  of  a  otep  input  and 
(2)  the  possibility  of  determining  the  moment  of  inertia  of 
the  aircraft  from  the  equations  of  notion. 

Transient  response  data  for  an  ?-oOA  airplane  Trere  obtained 
from  Cornell  Aeronautical  Laboratory,  Inc.,  (fir.  1  and  2).  The 
flight  testing  was  performed  by  Cornell  Aeronautical  Laboratory, 
Inc.,  and  information  pertaining  to  the  instrumentation  and 
method  of  tests  can  be  obtained  in  Ref«  3» 

The  present  report  is  limited  to  an  analysis  of  tv/o  sets 
of  data  consisting  of  the  response  due  to  an  impulse  input  to 
the  elevator  and  the  response  due  to  a  step  deflection  of  the 
elevator.  The  stability  derivatives  obtained  are  compared  to 
theoretical  values  and  those  previously  determined  from  steady 
state  oscillation  data  by  Cornell  Aeronautical  Laboratory,  Inc., 
(Aef.  3) 

The  loach  Ihimber,  altitude,  and  eg  position  are . essentially 
the  sane  for  each  run. 


G0H7BNTOTS 

ax.  a 

The  wind  axes  or  stability  axes  are  used.  The  origin  0 
lies  at  the  center  of  gravity  of  the  airplane.  The  OX  and  OZ 
axes  lie  in  the  plane  of  symmetry  and  OY  is  perpendicular  to  • 
it.  'The  OX  aids  always  points  in  the  direction  of  notion  or 
into  the  relative  .rind. 

X  axis,  or  longitudinal  axis,  po^it5_ve  forward 

Y  axis,  or  transverse  axi3,  positive  along  the  right  wing 

Z  axis,  or  normal  axis,  positive  downward 

LnidAIl  DIJPLACXilXiTS 

A  linear  displacement  along  a  positive  reference  axis  is 
considered  to  be  positive, 

ANGULAR  DISFLAC  !  '  MT3  AND    HIS 


An  angular  displacement  or  moment  v?hich  is  clockvrise  when 
viewed  from  the  origin  looking  along  a  positive  reference  axis 
is  considered  to  be  pooi,v\ve. 

7XLQCTTX-, ;  AXD  ACCUSATIONS 

Velocities  and  accelerations,  either  linear  or  angular,  are 
considered  positive  in  the  sane  sense  as  the  corresponding  dis- 
placements • 


COIITrtOL  SU:l?A3^  DlffLEC'JIOIiu 

Positive  elevator  angle  is  associated  with  a  downward  movc- 
ment  of  the  elevator  trailing  edge* 

TAIL  LOAD 

Positive  tail  load  is  associated  with  a  dov/nward  load  on  the 
horizontal  tail. 

UO  TATIOLI 

The  symbols  used  in  this  report  are  essentially  those  employed 
by  the  National  Advisory  Committee  for  Aeronautics.  Jymbols  per- 
tinent to  this  report  are  listed  below  Tor  convenience: 

V0  true  airspeed,  ft.  per  sec. 

q  pitching  velocity,  radians  per  sec. 

Iy  moment  of  inertia  about  1   axis,  slug  ft." 

Ky  radius  of  gyration  about  Y  axis,  ft. 

m  mass  of  airplane,  slugs 

W  weight  of  airplane,  lbs. 

c      mean  aerodynamic  chord,  ft. 

2 
g      acceleration  of  gravity,  ft.  per  sec." 

p 
3  wing  area,  ft. 

t  time,  sec. 

a,  incremental  normal  acceleration,  g's 

c^  wing  angle  of  attack,  radians 

<f  elevator  deflection,  radians 


£  dovmwash  angle,  radians 

Q  angle  of  pitch ,  between  X  axis  and  the  horizontal, 

radians 

p  mass  density  of  air,  slugs  per  ft.'' 

^  airplane  density  factor, M  -   m  /g  S   c 

rf  aerodynamic  time  unit,  sec.,  T  ■  a  /g  S  V0 

d  differential  operator,  da  d  /  d  _L  o  T  (d  /  dt) 

1^  tail  length,  ft,  (c»g.  to  2!?#  tail  chord) 

C^  initial  air-plane  lift  coefficient,  C^  =  V,r  /  v  /?  5  V 

Cjjj     airplane  moment  coefficient,  G^  »  v  /  jaS^c >  c 

2      ? 
h      moment  of  inertia  parameter,  h  -  21^.  //t^  c 

CL     slope  of  the  lift  curve,  \       Z   dCL  /  doC 

Gm     static  longitudinal  stability  criterion,  Qjg   =  dG,^  /do( 

G^,     rate  of  change  of  pitching  moment  coefficient  with 

rate  of  change  of  angle  of  attack  with  re- 
spect to  t/r 

C^^   airplane  damping  in  pitch 

Cja     elevator  control  power ,  Cn   »  ^Gjq  /j£ 

dfe/dct  rate  of  change  of  downwash  vrith  respect  to  ec 

i^     angle  of  incidence  of  the  tail,  radians 

Cm     rate  of  change  of  moment  coefficient  v.loh  respect  to  i-^ 

H*     observed  tail  load,  lbs. 

1      moment  arm  of  observed  tail  load,  ft,  (c,g#  to  $0% 

tail  chord) 

W*     weight  of  tail  outboard  of  strain  gage  fittings 

p 
cH-t     tail  loaci  coefficient,  Cjj.  -  11^  / :  p  V0  S 


rn 


■> 


V,  Subscript,  1,  denotes  that  the  stability  deriva- 
tives ere  for  the  "tail-off"  airplane, 

■ae 


m 


d<£ 


J 


THEOIU  AIID  ANAKSia 

lb3t  of  the  applications  of  mathematical  analysis  to  the 
motion  of  an  airplane  have  been  for  the  purpose  of  determining 
the  airplane fs  response  to  a  disturbance.  This  study  considers 
the  inverse  of  the  above  problem,  that  is,  the  determination  of 
the  stability  derivatives  of  the  airplane  from  the  longitudinal 
equations  of  motion  vfhen  response  information  is  given. 

In  setting  up  the  longitudinal  equations  of  motion  the 
following  assumptions  were  made: 

(1)  The  airplane  remains  at  constant  airspeed. 

(2)  The  equations  are  linear  with  constant  coefficients, 

(3)  Measurements  of  the  forcing  function  and  responses 
are  subject  to  random  errors  only. 

Because  of  the  first  assumption,  the  longitudinal  equations 
of  motion  for  the  short  period  transient  response  may  be  written 
in  terms  of  the  lift  and  moment  equations  only. 


(1)  ^^  ^  <f 


oC  4-doC  ..  6.9   -  ■_,<■  ir_  B  0 


(2)  C^  oC-f  ^cloCtC^  d9   -  hd20  s-  G^f 

Since  incremental  normal  acceleration,  n^,  is  more  readily 
determined  in  flight  than  angle  of  attack,  oO  ,  tlie  following  re- 
lation can  bo  utilized  to  incorporate  n^   in  place  of  oC  . 


(3)  <=c  =      1 


°L 


OC 
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Also,  sinco  the  pitching  velocity  is  usually  neasured  by  the 
rate  gyro  in  terms  of  q  ■  d  0  /dt,  the  following  relation  nay 
be  used  for  d  0  • 

(U)  d  9  m  Y  q 

Substituting  Equations  (3)  and  (1;)  into  Hqpxatlon  (l),  the 
lift  eauation  becomes 


^)     CL D8H-0LdiiatrOLocq-ll<fl. d<f  =     0 

2  It  /c 


ri 


ith  rig,  dn2,  q,  and  d  6   available  from  the  transient  response 


data.  (>r    can  now  be  determined.  CL,    can  also  be  determined 
fron  the  above  equation,  but  its  accuracy  will  be  limited  by  the 
fact  that  it  appears  in  a  constant  coefficient  that  is  quite  a 
snail  number. 

For  convenience,  Equation  (5)  is  regrouped  in  the  following 
form: 

(6)   ^  y  +■  diw  -  Cn<*     2. d<f  =  0 

where     (6a)  v  =    ^   iu  ■+  f  q 

2 

Because  of  the  relatively  poor  rate  gyro  measurement  of  q 
on  the  oscillograph,  as  shown  in  Fig«  1  and  2,  the  use  of  q  in 
the  solution  of  the  moment  derivatives  is  avoided.  By  solving 
for  d  9     in  Equation  (l)  and  substituting  d  B   and  also  oC  from 
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Equation  (3)  into  Equation  (2)  the  moment  equation  boco 

cL  0^ 

(7)      (CHcc+  -11  0mae  )n2  +  (Cndo6r  C^e  -  h  —l)dna 

C^  .         1  14 


CI  ^m  h 

1+  /c     CL         ^  Xt  /c     CL 


/c     UL 

In  order  to  avoid  the  inherent  inaccuracy  of  determining  second 
derivatives  from  the  transient  response  data,  Equation  (7)  is 
integrated  once.     The  resulting  moment  equation  is  as  follo-,vs: 


c  c  c 

(8)      C<^  + JL^*  )    I  r^d  Vr+  (C^*  C^p    -  h  J^-K 


-  hdnz  - f^_ _L  c^^  J:  GLcc)  J>d  t/r 


With  the  inertia  parameter,  h,  and  all  the  response  data  ioior.ii, 
this  equation  can  be  solved  for  the  "lumped  constants"  or  com- 
binations of  moment  derivatives. 

In  order  to  determine  uniquely  the  values  of  the  moment 
derivatives,  the  additional  equations  required  are  found  by  use 
of  the  observed  tail  load,  as  suggested  by  Itef«  k»     "Suppose 
that  the  horizontal  tail  of  the  airplane  were  totally  ineffective 
in  producing  any  aerodynamic  forces  or  moments  on  the  airplane. 
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Suppose  then  that  the  airplane  in  this  condition  were  subjected 
somehow  to  loads  at  the  tail  exactly  equal  to  the  tall  loads  ob- 
served in  any  practical  maneuver.  It  is  evident  that  the  response 
■would  be  exactly  the  observed  one,"  The  longitudinal  equations 
of  notion  for  the  airpla.te  in  this  condition  are 


(O)   ^QC 


■  +  doc     -  d&     -  J^  =0 


1 
(10)  ^m   +  Cq.    doLtC^,   d*-  hd^  =  -  —  Cnt 

1       1  1 

The  equation  of  ct  In  terms  of  n^  nou  becones 

/  i   -,     CL       \    Cnt 

(11    oC    =  — /  -^  -  nz 

CLoC      I      °L 

3ubstitutinrj  Equations  (!;.)  and  (11)  into  Equation  (?) 9   the 
"tail-off"  lift  equation  becomes 


(12)  ~ dynrq  +  iCn.-  %-     y 

CLc<  % 


1 
vrhere     (12a)   y  =  -—  C^  -  r^ 


With  all  variables  known  from  the  transient  response  data, 

C-l    ,  may  be  uniquely  determined. 
1 

By  solving  for  d  &  in  Equation  (">)  and  substituting  d  9     and 
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also  oc   .iron  Equation  (ll)  into  Equation  (10)  the  "tail-off" 
moment  equation  becomes 


Lc*. 


(13)  (C^  f  L  C^   )  7 


+  ((W  +  Si*  - h  -^- )  ^ 

1        1       2 


Si*    GL« 

nt 


-  hd2y  -  t  L  L  )  c 


2       °L 
h  GL«rn   i 


Since  the  derivative  Ora^,   is  small  enough  to  be  considered 

negligible,  it  is  omitted  hereafter. 

Once  again,  to  avoid  determining  second  derivatives  from  the 
transient  response  data,  Equation  (13)  is  integrated.  The  re- 
sulting equation  is  as  folio.  ,s: 


G-W     (Gmo,  +     ■  °<1    Si*    )  /  vd  t/r 


-+-  (Sie    -  h    Lo(l  )  y  -  **  = 

1  2 

°Lc(.  Si* 

i.  (    L    -JL)    /  cw  d'Vr 


CL  2  c 


yiH 


1U 


-where  w  -  (dy  -      .,.,---  Cn.    )   is  employed  for  simplification 

2GL 

of  calculations* 

In  the  laot  coefficient  of  Equation  (ll;)  the  term      1 

2 

is  snail  compared  to  l/c  and  is  omitted.  This  allows  for  the 

solution  of  h  and  hence  the  moment  of  inertia,  ly,  from  the 

flight  test  data. 

Jith  all  the  response  data  variables  known*  E     n   (lii) 

can  be  solved  for  h,  C^    ,  and  Gr ,  &      . 

1  1 

It  is  now  possible  to  set  down  six  equations  from  which 
the  unique  moment  derivatives  may  be  found. 


«J   ^-W-H^  t1"^. 


<16>  ^  aG-d^  4tS 

i 


(«)  ^  =  -^-4-  h  ~r 


GL<* 
(10)  G,^  -f-  Pels*  -  b     _  »  Constant,  coefficient 

from  Equation  (3) 


CLc 
n  .    -f-  ***      ^njg81  Constant,  coefficient  from 


do  a.    +    L*- 


2 
Equation  (8) 


1$ 


n 
(20)  l  a\         1       (Cn   .  ^   Ct  )  «  Constant,  coef- 
xt  /c   GL     ^   c      <*■ 

ficient  from  Equation  (0) 

These  six  simultaneous  equations  can  be  solved  Tor  the  six  ui>« 

knowns,  Cn^,  Cj^  ,  Gy^Q    ,  ^       ,  Cig.  ,  and  d£/doC  . 

The  solutions  of  the  equations  to  determine  Cl  ,  h,  the 
"tail-off"  coefficients,  and  the  combined  moment  derivatives 
all  depend  upon  the  use  of  redundant  data.  The  method  of  least 

squares  is  employed  as  the  procedure  in  handling  the  redundant 
data.  The  authors  are  well  aware  that  the  method  of  least  squares 
is  not  applied  in  a  rigorous  fashion.  Although  all  errors  in 
the  data  were  considered  random,  it  is  quite  possible  that 
systematic  errors  are  present. 

In  some  instances  variables  are  combined  to  form  new  va^- 
riables  (Equations  6a,  12a) •  Also  the  errors  are  assumed  to  be 
in  one  variable  only.  (Goc  Appendix  C.)  However,  it  is  believed 
that  application  of  least  squares  as  used  throughout  the  analysis 
is  a  reasonable  way  of  determining  the  stability  derivatives 
from  the  redundant  data. 
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RESULTS  AIID  DISCUSSION 

The  results  of  this  investigation  are  presented  in  the  tables 

below.  Results  of  previous  invest igations  from  both  theoretical 

analysis  and  flight  test  data  are  also  presented  for  comparison. 

^Analytical  Previous    Step  Function  ",K*Impulse 

Results    Flight  Test  Results  Function 

ilo suits     (Run  £301)  Results 

(Run  £063) 

i7CL              6.U3       6.19       6.22  £.91 

Cr               6.12       £.8        £.93  6.37 
*1 

h                                         0.0133              0.0l£0  0.0100 

Cjj^                                       0                    wide  scatter  -  0.0102  0.0?,76 
^^1                      '                                 -0.0121;  to 

4-0.2UO 

Gn^                                    0.££                  0.U8                  0.U3£  0.397 

^tf^Afce                   —               ~  °*Lt°3            -  0^39  -  0.1;6l 

GLcr 

Sloe"  Si©"*1                      ~  °-°90              -  0.1011  -  0.103 

i-(Cn  +  CL    h) -^0-°                 -Ul*16  -[^2^ 

It  /c  °L 

Gm  r                              -  0.63           -  0.36             -  0.96  -  0.86 

6 

^mcc                                -  °-2?            ~  °»35              -  0.321;  -  0.3£7 

C^^             -  0.012  '  -  0.022     -  0.023  -  0.033£ 

C^e               -0.021;     -0.032      -  0.014;  -0.0268 

Cnit             -  1.7U     -  2.U;      -  2.14  -  3.08 

d£/d«-            +  0.£l     +  0.63      +  0.683  +  0.7£6 


e.g.  ■  27;^  IlaC        U  ■  0.707 
■it- 
Approximate  averages  obtained  from  lief.  3. 

The  lift  due  to  the  deflection  of  the  elevator  was  omitted  in  de- 
termining moment  derivatives. 
M 
v   Derivatives  are  per  radian. 
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As  noted  iron  the  above  tabic,  the  results  of  this  investi- 
gation«  for  both  the  stop  function  and  the  impulse  function,  ore 

numerically  comparable  to  those  obtained  previously.  ?tobi  tliese 
results  it  appears  that  the  method  used  to  analyze  transient 
response  data  is  practical  and,  with  certain  precautionary  meas- 
ures to  be  discussed  throughout  this  section,  the  derivatives 
obtained  are  comparable  in  accuracy  to  those  determined  from 
frequency  response  methods. 

The  oscillograph  records  of  the  transient  response  data 
are  shown  in  Pic*  1  and  2.  All  values  of  the  basic  variables 

•oin  were  taken  within  a  one-second  interval  after  the  elcvato: 
deflection  was  applied.  The  one-second  interval  was  divided  into 
twenty  equally  spaced  points  at  which  the  readings  of  the  vari- 
able.-) were  taken*  The  readings  of  pitching  velocity  and  tail 
load  were  corrected  for  calibration  error  and  the  effect  of  tail 
weight  in  Appendices  A  and  B,  respectively. 

Che  rate  gyro  measurement  of  q  was  exceptionally  poor.  The 
oscil  atory  q  response  was  faired  as  shown  in  Fig*  1  and  2.  For 
this  i  eason  the  q  data  was  used  only  where  necessary v  that  is, 

in  the  calculation  of  C^   end  Cj.  . 

^  1 

-Iciypor.i  o  to  the  ^teo  -reflection  of  the  Elevator 


■    ihe  calculation  of  GT        «  GT        ,  h.  and  the  "tail-off" 

^->  oc  ^  or 

derivatives,  data  at  Points  3  through  20  were  used  in  the  least 

.  solution.  The  values  of  Ct   a  6.22  and  Cj,^  »  ^.03 

1 
agree  q  !.te  well  with  the  values  CL   -  6.19  and  C^  =5.0 

1 
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obtained  previously.  However,  the  inertia  parameter,  h  =  .31  , 
is  approximately  12  per  cent  higher  than  h  ■  i.  1133  as  calculated 
from  the  previously  lenown  value  oi  moment  of  inertia.  There  Ls 
no  apparent  reason  for  t :.e  discrepancy,  and  from  a  theoretical 
viewpoint  it  should  bo  possible  bo  obtain  sufficiently  accurate 
values  of  -he  Inertia  parameter. 

idix  D  contains  a  more  simple  equation  from  ..  dc  -  h 
o   bo  determined.  This  e  u;  bion  could  not  be  used  to  adv 
in  the      it  analysis  because  oi  its     -doncy  on  accurate 
raic  gyro  meaai  r      )f  pitching  velocity. 

Tne   combine      it  leriv;  bives  were  obtained  by  '^lio   solu- 
tion of  Equation  (G)  of  the  hi  lysis  section  of  this 
report.  Values  of  1   h-  sic  variables  were  tahen  at  the  even 
numbered  points  in  hi-;.  1.  The  nunbt  ?    >f  points  used  in  the 
least  squares  solution  '  as  reduced  bo  ten  In  order  bo  reduce  the 
>unt  of  calculations  involved. 

.hi  ins  cction  of  hi-;.  1  in  Leal  s       ie  slope  oh  the 
elevator  response  is  practically  nep  '  ;ible  subsequent  to  Point  2. 
A  first  attenpt  at  reading  the  slope  after   '  b  2  gave  a  maxi- 
mum and  minimum  d<f  =  £0. 000221  per  radian.      these  very  small 
values  of  d6   were  used,  completely  absurd  results  f  r  the  com- 
bined mo   b  derivatives  were  obtained.  However,  when  oho  value 
of  d  6     =0  was  assumed,  the  values  oi  the  combined  moment  deriva- 
tives were  reasonable,  hie  discrepancy  in  the  results  can  be 
the  fact  that  d  6     is  such  a  small  number  that  a 
very  slight  ;:.^:o^   in  lagnitude  may  bo  large  percentage  wise.  In 
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the  solution  or  the  simultaneous  equations  the  d  <S     tern  must 
be  eliminated  in  order  to  obtain  the  combined  derivatives.  As 
this  proces's  is  dependent  upon  either  a  division  or  a  multiplica- 
tion, it  is  them percentage  error  that  carries  through.  Also,  as 
the  entire  solution  depends  upon  the  subtraction  of  relatively 
large  numbers  to  obtain  a  suall  number,  a  large  percentage  error 
in  d«T   causes  completely  erroneous  results.  Therefore  the  as- 
>tion  of  d<f   =0  -was  used  to  calculate  the  tabulated  results • 
The  individual  stability  derivative-  are  in  close  agreement 
vs  h  those  previously  obtained  although  both  C^  .  and  C™    are 
so  r.-ihat  high.  In  the  calculation  of  the  moment  derivatives, 
the  exception  of  CQ  ,  the  value  of  C^,  g      tras  found  to  be 
rtant.  (This  applies  especially  to  the  impulse  data  and  is 
bro  ;ht  out  later  1X2  this  paper.)  The  previous  results  obtained 

fro:  oscillation  data  indicate  that  C^g   is  very  close  to  zero. 

1 
In  o  '  ler  "bo  determine  the  effect  the  magnitude  of  Cnvi-a   has  on 

1 
the  •  lues  of  the  moment  derivatives,  &tvq      was  assumed  to 

equal  -,ero.  The  table  below  indicates  those  results. 

°rad0   ■  -0.0102  Grade   ■  0 


0.   ,  per  rad.  -0.962  -  .963 

Vi    t  •       -Q.32J*  -0.339 

C-a  fc  -0.023  -0.028 

Gux^  -O.OUii  -0.039 

Gmit  -2.10;  -2.32 

d£/c  cc  0.633  0.719 
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In  general  the  results  obtained  from  the  step  function 
data  are  considered  to  be  good  representative  values  of  the 
stability  derivatives  of  the  F-GOA  aircraft* 

on^e  to  the  Ir.nulse  c flection  of  th-  -  leva tor 

The  values  of  Cj^  =  £.91  and  G^    a  6.37  do  not  agree 

with  either  previous  results  or  theory.  These  results  were 

obtained  from  an  analysis  of  Points  2  through  20,  excepting 

Point  10.  As  the  difference  between  Ct^.  and  Griar     is  due 

to  the  tail  lift,  it  vrould  be  expected  that  Gjr   would  be  greater 

than  C-r  „    •  A  possible  explanation  might  be  that  the  poor  pitch- 
*1 

ing  velocity  data  caused  the  discrepancy.  As  the  difference  be- 
tween Cl   and  Cj,    should  be  very  snail,  an  average  value  of 


1 
^har   =  ^Lar  a  6«13  Tras  assumed  for  use  in  subsequent  calcula- 
tions. 


The  values  Df  h  and  the  "tail-off"  derivatives  -..-ere  found 
to  vary  somewhat  from  the  results  obtained  by  Gorncll  Aeronautical 
Laboratory,  Inc.  In  particular  h  «  0.J10  was  obtained  as  com- 
pared to  h  =  0.01.33  obtained  previously,  and  h  =  0.015>  as  deter- 
mined from  analysis  of  the  step  data,  ^ac    f   --n  this  case, 

changed  sign  to  C^jg   »  0.0176.  This  value,  although  well 

1 
within  the  range  of  the  scatter  obtained  by  frequency  response 

analyses,  causes  some  trouble  in  the  determination  of  the  actual 

derivatives.  This  point  is  discussed  later  in  this  section  at 

somewhat  greater  length. 

In  the  case  of  the  impulse,  as  well  as  with  the  step,  snail 
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tcms  were  troublesome.  A  first  attempt  go  analyze  the  data  at 
points  evenly  divided  throughout  the  range  (even  numbered  Points 
2  through  20)  ended  in  absurd  answers  for  the  combined  constants. 
However,  by  making  the  assumption  that  the  lift  due  to  the  de- 
flection of  the  elevator  as  negligible,  ihe  tabulated  results 
were  obtained.  The  above  assumption  is  reasonable  because  the 
elevator  is  deflected  for  a  very  short  time  only.  The  us,e  of 

the  assumption  in  effect  causes  d<f  3  <f   »  0  for  all  points, 

« 
and  thus  eliminates  the  diffic  lties  of  snail  numbers,  as  pre- 
viously discussed. 

The  effect  of  Cnd0   on  ^-ie  calculations  of  the  stability 
parameters  is  quite  pronounced.  If,  as  v„ras  done  in  the  consid- 
eration of  the  step  function,  Gm^   -  :'-^  assumed  to  be  G-  y  *      =  0, 
the  derivatives  change  in  the  manner  indicated  below: 

Previous     Grid^   a  0     0^^   a  0.0176 
Results         l  i 

Cn  ,  oer  rad.  —0.36  -  .36  -0.36 

Cn^  ,  per  rad.  -0.33'  -  .331  -0.3^7 

C^^  ,  per  rad.  -0.022  -0.02$  --.0335 

C^e  ,  per  rad.  -0.032  -.03^1;  -0.0263 


CT1       ,  per  rad.  -2.  h  -2.1*6  -3.:jG 

1t 

dt /doc  0.63  0.701;  0.7>6 
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The  above  table  indicates  that  possibly  the  value  of  ^ix^ 

1 
should  be  closer  to  zero.  At  any  rate  it  is  evident  that  the 

success  in  the  method  of  analysis  used  is  dependent  upon  the  ac- 
curacy of  Cj^   . 
1 
A  comparison  of  the  results  obtained  from  the  transient 

data  with  those  previously  obtained  from  frequency  response 
data  indicate j  that  possibly  the  step  function  data  is  best 
suited  for  this  type  of  analysis.  This  result  is  as  expected, 
as  the  equations  used  depended  to  a  great  extent  upon  the  values 
of  integrals.  The  step  function  data  yielded  larger  integrals, 
and  the  errors  in  measurement  are  smaller  percentage  ".rise.  How- 
ever, if  the  equations  of  Appendix  D  are  utilized,  the  impulse 
nay  be  as  good  as  the  step  function. 

/mother  consideration  might  be  that  the  frequency  content 
of  a  step  function  is  quite  low  as  compared  to  that  of  an  impulse 
in  the  high  frequency  range.  This  has  the  advantage  that  any 
higher  order  effects  (important  at  high  frequencies)  that  were 
neglected  from  the  equations  of  motion  vrill  be  unimportant .  How- 
ever, the  frequency  content  of  a  step  function  is  very  high  in 
the  low  frequency  range.  As  phugoid  frequencj.es  are  also  Low, 
precautions  must  be  observed  to  sec  that  the  assumption  of  con- 
stant airspeed  is  valid. 
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GozicLusioiTj  a::d  mQomsmkTTO] 

The  results  of  the  present  investigation  lead  to  the  fol- 
lovfing  cone  liis  ions: 

1«  It  is  feasible  to  obtain  the  longitudinal  stabil- 
ity derivatives  of  an  airplane  by  the  "Equations  of  Ibtion  Ifethod 
of  Analysis  of  Transient  Response  Data," 

2.  The  step  function  data  yields  better  results  than 
impulse  data  if  the  equations  used  depend  upon  integrals,  as 
they  did  in  this  study.  This  is  not  necessarily  true  if  varia- 
tions of  the  equations  of  motion  ere  used  in  ■which  integrals 
are  not  involved, 

3,  With  good  instrumentation  it  should  be  possible 

to  obtain  the  inertia  parameter  h,  and  thus  the  moment  of  inertia 
of  the  airplane  by  transient  response  analysis, 

lu  Reasonable  values  of  the  longitudinal  derivatives 
of  the  F-80A  for  a  liach  IJumber  of  H  «  0.707  and  a  eg  position 
of  27^  ma   are: 

Cl  6,22  /  radian 


C-  -0.96  /  radian 


ra  oc 


-0.3214  /  radian 
C^  -0.023  /  radian 

Coa  q.  -O.Olilj.  /  radian 

J?«  The  use  of  the  BEquations  of  Ibtion  Ijethod"  is 
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dependent  upon  small  differences  of  large  numbers  so  that  va- 
riables of  very  small  magnitude  ssu*t  bo  avoided  v/here  possible. 

6.  The  accurate  evaluation  of  Cjj,    is  very  important 

1 
in  the  subsequent  calculation  of  the  moment  derivatives  in  the 

method  used  5_n  this  investigation, 
ilecommenda  t  ions : 

1.  Additional  flight  data  should  be  analyzed  by  J.ie 
method  outlined  in  this  reoort  in  order  oo  develop  the  method 
further.  This  additional  data  should  include  information  at 
various  !Tacn  llumbers  and  center  of  gravity  positions. 

2.  The  equations  outlined  in  Appendix  D  should  be 
utilized,  provided  the  instrumentation  gives  sufficiently  ac- 
curate pitching  velocity  response. 
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TABLE  I 


DBE2JSI0    ;  AI1D  LEADING  PARTICULARS 

P-    )A  Serial  No.  1^169 


I.  Principal  Di-cnuions 

A.  Airplane  -  Gone    al 

Span 

Length  (overa]  ) 

Height 
B#  .Ving 

Aspect  Ratio 

Lean  Aerodyna  .c  Chord  (MAC) 
C  •  Empennage 

Tail  length,  ti 

Tail  length,  3  | 

Aspect  ratio 

?x>ot  chord 

Tip  chord 

Elevator  hinge  center  line 

II.  Areas 

rtng 

Total  v.-ing  area 

B.  2-rpennage 

Total  horizontal  tail  area 


Di  tensions 

38 «  10i  in. 

3ii»  6  in. 
11 «  k  in. 

6.38 
0O.6  in. 


lk.87  ft. 

15.338  ft. 

5o6 

li.33  ft. 

1.583  ft. 

75>/j  chord 

237.6  sq.  ft. 


U3«5  sq»  ft. 
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'LJSL-1   I 
(continued) 


Dimens  ions 


III.  Uass  Characteristics 

A.  Yfeight  o.t  aircraft 

Ruii  ^301  10,51$  lb. 

.-1  5068  -  10,263  lb. 

honent  of  inertia  about  '1     axis 
Run  $301  l.:,130  slur  ft.2 

Run  ^068  1^,130  slu-  ft.2 

C.  C.G.  Position,  %  ' 

Run  5301  27.' 

Run  5063  27^ 

D#  V/eight  of  tail  outboard  oT  strain 

garo  installation,  w-fc  100  lb. 

IV,  Oscillograph  Trace  Sensitivities 
wane  for  both  5063  and  5301) 

&  0.i£5  dec.  /  in. 

q  Ii.39  deg«  /  sec.  /  in. 

ns  0.351i  3  /  in. 

llt  151  lb.  /  in. 

V.  Other  Particulars 

.'-.  Run  53  'lj  Step  function  data 

V0  726  ft.  /  sec. 

p  0.001271  slugs  /  ft. 3 

Iiach  ;?a:nber  0.1 

/atitude  20,lS0  ft. 
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TAKLE  I 

(continued) 


Dimensions 

r 

l»ii9  sec. 

M 

160.8 

CL 

0.132 

B. 

Run  5o60,  Impulse  function  data 

V0 

733.5  ft.  /  sec. 

e 

0.0012>3  slugs  /  ft.3 

Mach  1  lumber 

0.703 

Altitude 

20,230  ft. 

r 

l.U^2  sec. 

/>c 

158.6 

cL 

0.1273 
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1,  3tep  Function  Response  Data,  Run  £301 30 

2.  Impulse  Function  Response  Data,  Run  5>068  ...  .  .31 

3«  Variation  of  y  with  Tirae^  Run  £301 32 

lw  Variation  of  y  with  Time,  .iim  Jo6o 33 

£•  Variation  of  Tail  Load  Coefficient  with 

rime,  Run  5301 3h 

6.  Variation  of  Tail  Load  Goof  fie  ".out  with 

WiBB,  :cun  5063 3^ 

A]  roendix  A 

1A«  Rate  Gyro  Frequency  Response 39 
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Constant •  U0 

3A,     Variation  of     J qdt  with  rime,  Run  >3--l  •  •  •  •  ill 
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APPENDIX  A 


A  METHOD  FOR  GQHVEROTG  RATS  GIRO  FBEQUE8CE  BESFQHS8 

CALIBRATION  DATA  TO  A  FORM  APPLICABLE  TO  TRAIIJISIIT  INPUTS 


Theory 

Because  of  the  fact  that  the  rate  gyro  used  to  measure  the 
pitching  velocity  was  calibrated  vrith  a  sine  wave  input,  it  was 
necessary  to  convert  the  calibration  data  in  such  a  manner  that 
it  could  be  used  for  transient  inputs.  The  frequency  response 
calibration  data  are  presented  in  Fig.  1A  as  amplitude  rate  and 
phase  angle  versus  frequency  curves. 

It  is  assumed  that  the  rate  gyro  is  essentially  a  second 
order  system  and  can  be  defined  by  the  equation 

(1)  qts^r'^  ir^+% 

or  (2)       a  -  G  (X+rrM-  X      -^g-    *  ) 

%  wn 

where 

wn  s  natural  frequency  of  the  system 

?  ss  the  damping  ratio 

%  "  pitching  velocity  indicated  by  the  rate  gyro, 
radians  /sec. 

q  *  pitching  velocity  of  airplane,  radians  /  sec. 

C  m   sensitivity  of  the  rate  gyro,  radians  /  sec.  /  in. 

X  =  displacement  of  oscillograph  record,  in. 
Therefore  if  /£  and  wn  for  the  system  are  known,  it  is  possible 
to  determine  the  pitching  velocity  of  the  airplane.  In  the  present 
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analysis  Equation  (2)  was  integrated  and  the  integral  of  q  was 
obtained.  This  was  clone  in  order  to  avoid  talcing  second  slopes 
01  idle  recorded  data. 

(3)    /qdt  s  -V-  <L  -f  -^  n  +   fc   dt 

If' Equation  (l)  is  transposed  to  operator  Tom  and  solved 
Tor  the  amplitude  ratio  of  qQ  /  q  the  following  equation  is  ob- 


Then  by  applying  the  well  known  "jit"  substitution  for  the  operator 
"s"  the  equation  becomes j 


(5)  AR(jw)  -  ..2 

(l-j£.)*<2fJ~  )j 

>Tn         n 


and      (6)  |  AH  1  ■ 


f(l-^-)2^(2fJI-  ): 


The  phase  angle,  d   ,  is  determined  as: 
(7)  <j>    =  -  arc  ton  2  ^  ^rT" 

i-(f-)2 

n 
Tlie  dancing  ratio,    £  ,  for  the  system  can  be  obtained  by 

comparing  the  plot  of  amplitude  ratio,  Ad,  versus  phase  angle,  <j>  , 
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to  constant  damping  ratio  curves  of  a  standard  second  order 
system  (Fig*  2A).  If  constant  damping  ratio  curves  are  not 
readily  available,  they  can  easily  be  plotted  from  Equations 
6  and  7  by  assuming  values  of  w/wn  and  solving  for  the  corres- 
ponding values  of  AR  and  <f>   at  the  desired  damping  ratio. 

:-'ig.  2A  indicates  that  the  damping  patio  for  the  present 
system  is  approximately  -%  »  0«l£.  Because  of  the  wide  scatter 
of  points  it  is  difficult  to  fair  the  data.  However,  as  the 
pitching  velocity  correction  is  small,  a  small  error  in  $\,  does 
not  materially  affect  the  resulting  pitching  velocity  of  -the 
airplane. 

As  £  is  now  hnown,  it  is  possible  to  assume  a  value  of  wn 
and  to  determine  both  Ail  and  S   for  each  of  several  frequencies 
(Equation  6) . 

A  comparison  of  the  values  of  AH  and  <f>  to  points  on  the 
curve  in  Fig.  1A  indicates  whether  the  assumed  value  of  wn  is* 
the  correct  one.  In  the  present  case  wn  ■»  21.2  fits  the  curve 
quite  vrall  (Fig.  2A) . 

If  the  values  of  £  =  0.l£,  wn  -  21.2  radians  /  sec,  and 
C  ■  0.0767  radians  /  sec.  /  in.  are  substituted  into  Equation  (2) 
and  the  equation  integrated,  the  following  expression  for  the 
correction  of  the  rate  gyro  reading  is  obtained* 

(G)  /qdt  B  0.0767  /xdt  V-0*00108£X  +  0.0001702X 
A  plot  of  /qdt  versus  time  is  shown  in  l?ig.  3A.  The  final 
values  of  q  were  found  by  reading  the  slopes  of  this  curve  at  the 
points  desired  in  the  calculations. 
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APPiSIDIX  B 


horizontal  tail  load  correction 

FO:i  MEETIA  LOADS 


The  tail  load  response  data  shown  in  Irig.  1  and  2  were 
measured  by  a  strain  gage*  The  response  shown  is  composed  of 
both  aerodynamic  loads  and  the  inertia  loads  of  the  stabilizer 
outboard  of  the  strain  gage  installation*  The  measured  tail 
load,  (llt)m>  H^st  ^e  corrected  as  follow: 


(1)  Ut  =  (IIt)n  -  «,  £  ^  -  Jfl_] 


where  w**  «  weight  of  tail  outboard  of  strain  gage  fittings. 

For  the  F-30A  strain  gage  installation,  ify  «■  100  lbs,  and 
1  =  15.332  ft.  Hence  Equation  (l)  becomes 

wts  (1Tt)m+100nz+U7.6  q       . 
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appendix  c 

SAMPLE  CALCULATIONS 

A  complete  outline  of  the  calculations  for  the  stability 

derivatives  for  the  case  of  the  elevator  step-function  input 
(pjun  3301)  is  presented  below. 

The  values  of  the  basic  variables,  for  any  one  instant, 
were  taken  from  Pig*  1  and  are  listed  in  fables  I  3  -  IV  E. 
All  measurements  of  the  basic  variables  were  determined  within 
a  period  of  one  second  after  the  elevator  deflection  was  applied. 
The  variable  q  was  corrected  for  instrument  calibration  error 
as  outlined  in  Appendix  A.  The  variable  :!-£,  was  corrected  for 
the  presence  of  horizontal  tail  inertia  loads  as  outlined  in 
Appendix  B, 

The  slope  of  the  lift  curve,  Cj,   ,  was  determined  as  folio. s: 


-   v  y-  diL;,  -    6       dS    =  0 


cL        "it '/ 


where 


v 


nz  -f  rq 


Applying  the  method  of  least  squares  and  assuming  that  the 
greatest  error  lies  in  v,  the  normal  equations  are 


^L.  y  (cLnz)2  -  ^S y  cU  .  d^-?-^  .  dna  -  0 


CL    K~,r    .      Or. 


Y*£    •  dnz  -  *XS  Y~(d«T  )2J-Yy   .  di"  =  0 

^- —  !•:•,  /c  Ct   /- —        ^- — 


T-<*  *~  H  /c  cL< 


The  numerical  values  for  the  summations  indicated  from 

Tabic  V  2  can  now  be  substituted  and  the  equations  solved  for 


°L  /  °iv  ' 


CT  /  CT    a  0.021223 


si 


incc     Cj,  «  0.132 


then     G^   =  6.22  per  radian 


'Hie  slopo  of  the  lift  curve  for  the  "tail-off"  condition, 
Ct    ,  was  determined  as  follows: 


£l_   dy  =  -Tq^J  Cnt  -   ^   y 


^1 


1 

where    y  =  (  -*--  Gn  -  nz  ) 

A  graph  of  the  y  variable  versus  time  is  show:  :_:.  3*  This 
graph  was  used  for  determining  the  values  of  dy  for  any  one  in- 
stant (or  point). 

Applying  the  method  of  least  squares  and  assuming  that  the 
greatest  error  lies  in  7"q>  the  normal  equation  becomes 


£dy  -  rci  -  -|— -£  W)V£<fr  •  CS<W  -J§r  •  (%-y)  = 


The  numerical  values  for  the  summations  indicated  from  Table 
71  E  were  substituted  and  the  equation  solved  for  0^  /  G-,^  . 


k6 


GT  /  Ct    =  .022261 

1 


CL^  =  £»?3  Per  radian 


The  values  01  h,  Cn    ,  and.  C™,^  vrere  determined  as 

1        "  1 
follows : 


$*«+■   ^"1   °>M»  )  .yyc!  t/^C^  -  h  j^l) 


cnt  d  t/r 


-  hw  fe»                1 

1 

cL 

c 

1 

7  =       CL       Oat  - 

nz 

where 

w  •  dy  -  Cl<1   Cnt 
2  oL 

Graphs  of  the  variables  y  and  C^  are  shown  in  Fig.  3  and  %, 
respectively.  These  graphs  vfere  used  in  dote.     1 ;  the  integrals 

of  the  variables* 

Putting  the  above  equation  in  the  form 

A  /yd  t/r/By  -  hw  =  -  C  Jcnt  d  t/r 

and  applying  the  method  of  least  squares $  assuming  that  the  great- 
est error  lies  in  v,-,  the  normal  equations  become: 

X  v,r  •  /°nt d  ^^ZJ°»t  d  ^ )2  -  irz^d  ^  • 

yOnt  d  t/r   -  -^rjj  •  /Cn,.  d  t/r  =  0 


hi 


-  *j-^y  •  f yd  t/r  *  o 


f!"2-' 


The  numerical  values  for  the  sumat  ions  indicated  iron 
Tables  VII  E  -  T.Z  S  were  substituted  in  the  equations  and  the 
equations  solved  for 


C/h  ■  -  6,6Ljlu7 
A/h  =  30.4J4.66 
B/h  -  -  3.7299U 


".There 


CT       C       0,132 
therefore,  h  =  O.Ol^Ol. 

The  value  of  h  determined  from  ground  test  of  the  inertia  of  the 
airplane  was  found  to  be  h  »  0.0133  and  is  used  for  all  further 
calculations  • 

CL, .  _   and  CL,     are  solved  for  as  f ollov.-s : 


B  =  C^  -  h  L*±        z  -  3.72994  h 
1       2 


U3 


C^^      -  0.0133     (  ~\ )  =  -  3.72991;  (0.0133) 

1 

(V,  =  -  0.0102  T>er  radian 

1 


cL 

A  e  W  -f       "V    Pnd^    =  30.U466  h 
1         2  1 


Gm      +        n        (-  0.0102)  -  30.1U+66  (0.0133) 

1 

CL  -  0.1<3£l8  per  radian 

m<X. 


The  values  of  the  combined  moment  derivative  coefficients 
v/ere  determined  as  follows: 


_£ (Gn   V-GL      -~-  )    /f  d  t/r 


0  p 

_*   •  2_ (c,^ x <^  )(f ^ _^  *    <tf 

It  /c         CL  -H  /°         GL 


=  hdng; 


I?or  the  step  function  elevator  deflect  ion  the  derivative  d£    vras 
talcen  as  dS      «•  0  at  the  points  in  question. 
Putting  the  abovo  equation  in  the  form 

DJnz  dt/r-fE%*  F  fs  d  t/r   -  g£    a  hdnz 


U9 


and  applying    he  method  oJ  least  squares,  assuming  that  the 
greatest  error  lies    Ln  dr^,  tiie  normal  equations  become 

^  Mn,  .    Jnj-d  t/r      -  D  /_(  /nzd  t/r  )2-  e£x%  .  j^  d  t/r 

-j~  ?^  fid  Vr     .  ynad  t/r-rGy  f«  fa*  d  t/r  =  0 
^  **z  •  "-idnz  -  D<^nz  •  y"na  d  t/r  -  S^(nj|)2-f  ^/na  .JXd  t/r 
-f-    (jjng  .  d£   =  0 

j)  hdnz  •    /<f  d  t/r    -  D/  fnz  d  t/r  •  /Xd  t/r  -  E/ n^  .J<f&  t/r 

JTf •  hdn,  -  'D^f.fnz  d  t/r    -  sTf.  r^fF^.  /Jd  t/r 

-f    G^  £   )2  =  0 

After  substituting  the  numerical  value  for  the  su         Lons 
as  indicated  from  Cables  X  J]  -  XIV  E,   the  follovfing  constants 
rrcre  determined: 


D  •  -     .;/.l  -     Cm  +     ^  *      A 


E  S  -  0.108  «     C^  -h  C^       -  h       ^L<* 


SO 


^JT  1  Xt 


Five  of  the  dynamic  stability  derivatives  were  solved  from 
the  sinultaneous  equations  listed  below* 


^     =     ^     G*it     l1'^) 


1  It 

3nd<9-  "  ^*d<9  "^  ~^T    "IT"    ^it 
1 


1        It  d6 

^oC  "      "sZ       C  -t     doC 


CL 
%  + —  W    --0jj6l 


Subs  Lng  the  known  values  in  the  above  equations,  they  becone 


°v  =  °-^^cn.t  (i--4i-) 


C^   =     -  0.0102  -f-O.01301  Cn± 


Cm**    -    0.01381  0n±t  — ;- 


^n^K    -/-    G^lj9.      a  -  0.066 7 
Cn  -/-  3.11  Cn^     =  -  0.l|6l 
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Solution  of  tiie  simultaneous  equations  yields 

Cm         o  -  0#32i;2  Dor  radian 

C^    =  -  0.023  per  radian 
C^g^  =  -  Q»0hk  per  radian 


C--H.    =  -  2mhh  per  radian 
"  -t 


de/dod  =  0.683 


Tho  value  ox  C*-    v/as  calculated  from  the  solution  of  the 
combined  moment  derivative  coefficient,  2,   determined  previously 


Cm 


o  op 


Jl_  )  /  0.32U  *  6.22   (2.22)7   =  -  Ww2! 

0.132   ^  J 


C_   ■  -  O.962  per  radian 

S 

The  solution  01  the  stability  derivatives  in  the  case  oj?  the 
impulse-function  elevator  deflection  was  carried  out  in  a  similar 
manner. 
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jndix  d 


.  VARIATION  IN  TIE  FORM  OS     EQUATIONS  OF  LDTION 
FOR  DETERMINING  DINA?  IG  STABILITY  DERIVATIVES 


The   ie  of  the  variable  q  was  eliminated  in  this  investiga- 
tion where  ver  possible  because  of  the  poor  rate  gyro  measurement 
of  its  trr  sient  response.  However,  by  maximum  use  of  q,  the 
use  of  hig  r  derivatives  as  variables  in  the  moment  equat:'  ons 
can  be  avoi  3d  and  the  number  of  variables  reduced. 

Squati  1  (3)  in  the  THEORX  AND  ANALYSIS  section  of  this 
report  can  1   replaced  by  the  equation 

<V    Ji  7    -rCmde?   -rhdq--i-Ont  =  0 
%.x  1 

from  -which  h,  m   ,  and  C™    can  be  readily  calculated. 
Equation  (  .lj)  can  also  be  replaced  by  the  equation 

(C    _^L  -  °n«  A^  }«|  +  (  5l£-   ^  .  C  ,  )  <f 

***   t  cL0C       Vc  "^    J 

+  f*W+  C^^q-hrdq^O 
from  which  new  combined  moment  derivative  coefficients  can  be 
calculated. 

If  a  good  rate  gyro  is  used  to  obtain  the  flight  test  x^espon^e 
of  q,  these  equations  may  give  more  accurate  results  than  the 
equivalent  equations  that  ere  utilized  in  this  report. 


c! 
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APPENDIX    E 

Tx\BLS  IE 

b; 

toIC  VARIABLES 
Ran  5301 

Point 
Irumber 

Time 
(sec.) 

%<g) 

**% 

Tcf  (racl.) 

v" 

1 

0.05 

-0.01061 

-0.739 

-0.011921 

-0.0126 

2 

0.10 

-0.01U7 

+  1.056 

-0.0283 

-0.0292 

3 

0.15 

+0.03135 

1.679 

-O.OU025 

-0.0331 

It 

0.20 

0.0956 

2.376 

-0.0539 

-0.0526 

5 

0.25 

0.1877 

2.850 

-o.o585 

-0.0U611 

6 

0.30 

0.283 

2.350 

-0.0760 

-0.0573 

7 

0.35 

0.379 

2.850 

-0.0335 

-0.0535 

3 

o.Uo 

0.U71 

2.690 

-0.0879 

-0.0568 

9 

o.U5 

o.^ 

2.585 

-0.0909 

-0.05U0 

10 

0.50 

0.6U1 

2.265 

-0.095U 

-0.0531 

11 

0.55 

0.720 

2.110 

-0.095U 

-0.01*79 

12 

o.6o 

0.791 

1.951 

-0.095U 

-0. 01.32 

■     13 

o.65 

0.850 

1.635 

-0.0951; 

-0.0393 

111 

0.70 

0.9011 

UhZk 

-0.095U 

-0.0358 

15 

0.75 

0.9U6 

1.137 

-0.0909 

-0.028I1 

16 

0.30 

0.935 

I.003 

-0.031*9 

-0.0198 

17 

0.35 

1.012 

0.792 

-0.08U9 

-0.0181 

18 

0.90 

1.033 

0.739 

-0.03U9 

-0.0173 

19 

0.95 

1.060 

0.528 

-0.0835 

-0.0135 

20 

1.0 

1.080 

O.369 

-0.0320 

-0.0107 

<gt-    % 

+    T^Q 

• 

5U 


Point 

limber 


TABLE  II  E 

■IC  VAi:7 
am  5301 


Time 
(sec.) 


d6 


(0^2)7  Cxif/2 


1 

0.0$' 

-0.00163 

-0.0001076 

-0.000308 

2 

0.10 

-0.02368 

-0.001563 

-0.0021*95 

3 

0.15 

-0.06125 

-o.ooUolt 

-  .00191* 

b 

0.20 

«f  0.00221 

-0.11510 

-0.00761 

-0.001315 

5 

0.25 

0.00221 

-0.19930 

-0.01320 

-0.000793 

6 

0.30 

0.00221 

-0.2872k 

-0.01896 

-0.000230 

.  7 

0.35 

0.00221 

-0.37597 

-0.02U8 

+0.000200 

8 

O.IiO 

0.00221 

-0.1&&9 

-0.0306 

0.0001*95 

9 

o.U5 

0.00221 

-0.5l4.73b 

■^0.0361 

0.000302 

10 

o.5o 

0.00221 

-0.62517 

-0 .01*13 

0.0010U5 

11 

0.55 

0.00221 

-0.70105 

-0.01*63 

0.001250 

12 

o.6o 

0.00221 

-0.76955 

-0.0503 

0.00H415 

13 

0.65 

0.00221 

-0.82620 

-0.05U6 

0.001570 

Hi 

0.70 

0.00221 

-0.87770 

-0.0530 

0.001733 

15 

0.75 

0 

-0.91715 

-0.0606 

0.001903 

16 

0.60 

0 

-0.95U05 

-0.0630 

0.00201*5 

17 

O.'Jp 

-0.00221 

-0.9796 

-0.06U7 

0.002138 

10 

0.90 

-0.00221 

-1.00U7 

-0.0663 

0.002200 

!9 

0.95 

-0.00221 

-1.0258 

-0.0675 

0.002260 

20 

1.00 

-0.00221 

-1.0U58 

-0.0690 

0.002261 

\ 


-  n. 


?5 


III  E 

BASIC  VARIABLES 
ilun  5301 


?o5_nt       jMjbb  dy 

toer     (sec*-) 


Jyd 


r 


TT 


/Ontdt.  ^ 


dt 


3 

h 
5 
6 

7 
3 

9 
10 
11 
12 

13 
Hi 
1$ 

16 
17 
18 
19 
20 

w 


0.05 

0.10 
0.15 

0*20 
0.25 
0.30 
0.33' 
O.liO 

o.U5 
0.50 

o.# 

0.60 
0.65 
0.70 
0.75 
0.30 
O.05 
0.90 

0.95 
1.0 


-0.3125 

-0,760 

-1.31:1 
-1.378 

-2.655 
-2.6?5 
-2.655 
-2.655 
-2.I460 
-2.232 

-2.085 
-1.81.8 

-1.6h0 
-1.355 
-1.191 
-O.969 
-0.319 
-0.671 
-0.531 


-0.lw6 


-0.00002015 

-0.0005062 
-0.0019162 
-O.OOU066 
-0.010066 

-0.013126 

-0.029206 

-0.0U3296 

-O.O6O386 

-O.03O1U6 

-O.IO2366 

-0.12697 

-0.15372 

-0.13222 

-0. 212 37 

-0.2U37U 

-0.27621- 

-3.30961; 

-0.3U33U 

-0.3736U 


-1.3189 
-2.6192 
-2.61a6 
-2.66UO 
-2.6773 
-2.1;960 
-2.2790 
-2.1112 
-1.9116 
-1.7106 
-1.14328 
-1.2766 
-I.0610 
-0.9150 
-O.7698 
-0.6825 
-0.51:76 


-0.0000168 

-0.0001258 

-0.0001763 

-0.0002861 

-0.0003552 

-0.0003338 

-0.0003927 

-O.OOO3692 

-0.0003256 

-0.0002635 

-0.0001867 

-0.0000971 

0.0000051 

0.000117U 

0.OQO2l|O3 

0.0003722 

0.0005115 

0.0006567 

0.0003059 

0.0009739 


-0.0003150 

-0.0011215 

-0.002231l 

-0.003927 

-0.005996 

-0.00357U 

-0.01131U 

-0.01U250 

-0.017203 
-0.020210 

-0.023268 

-0.023363 
^).029la3 
-0.032liliQ 
-0.035U32 

-0.03333^ 
-0  .01^266 
-0.0iUil52 
-0.0U69lil 
-O.OU9630 


dy 


2CL 


°nt 


TABLE   IV  E 


% 


BASIC  VARIABLES 
rtun  5301 


Point         Time  ^(rad.)  J$d  t/r 

Number       (sec,) 


hdi^ 


J  njjd  t/f 


2 

0*10 

0.01002 

0.0003555 

o.oiia 

-0.000712 

it 

0.20 

0.01039 

0.001057 

0.0315 

-fO.OOll^O 

6 

0.30 

0.01039 

0.001759 

0.0379 

0.01U39 

8 

o.Uo 

0.01061 

0.002U60 

0.0358 

0.03981 

10 

O.pO 

0.01068 

0.00316^ 

0.0301 

0.07711 

12 

0.60 

0.01068 

0.003890 

0.0260 

0.12513. 

Hi 

0.70 

0.01103 

0.001*620 

0.0189 

0.13171 

16 

0.80 

o.oim 

0.005360 

0.0133 

0.21+511 

13 

0.90 

0.01090 

0.006100 

0.0098 

0.31231 

20 

1.00 

0.01090 

0.0063UO 

0.00U9 

0.33361 

PRODUCTS  AND  SUIB  OF  BASIC 

FOR  DET  LTCTATIOH  OP  Gt 

Run  5301 


57 


Point 

Purtber 

i 

vdS 

(d«)(dnz) 

(dn2 

(dnz)2 

h 

-0.12i.9pl; 

-0.00011622U 

0.00525096 

o.oooool;83ia 

5.6;6'38 

5 

-O.I3I4I4 

-O.OOOIOI903 

0. 00629850 

o.oooool;88la 

.12250 

6 

-0.163303' 

-0.000126633 

O.OO629350 

o.oooooiiosia 

G. 12250 

7 

-0.166723' 

-0.000129235 

O.OO62985O 

0.0000048841 

8.12250 

8 

-0.152792 

-0.000125523 

0.005914I90 

.  .'Oooouscip. 

7.23610 

9 

-O.139&0 

-O.0001193U0 

0.005712 35 

0 .oooooU33la 

6.63223 

10 

-0.120272 

-0.000117351 

0.00500565 

o.oooooli33ia 

5.13023 

11 

-O.IOIO69 

-0.000105359 

O.OO4663IO 

O.OOOOOl;831a 

k.h$210 

12 

-O.o0UlG3'7 

-0.00009x362 

O.OOl>31171 

o.oooooU38la 

3-O06I4O 

13 

-0.06U2 555 

-0.000036853 

0.00361335 

.OOOOOuSSla 

2.67323 

lU 

-0.05:09792 

-0.000079118 

0.003li'J0li 

O.OOOOOlt331a 

2.02773 

15 

-0.0337108 

0 

0 

0 

1.U0397 

16 

-0.019859U 

0 

0 

0 

1.J0601 

17 

-0.01U3352 

-fO.OOOOh.OOOl 

-0.00175032 

O.OOOOOli881a 

0.62726U 

18 

-0.012  7  8I-7 

0.000038233 

-0.00163319 

0.000004881a 

.  '46121 

19 

-0.00712oO 

O.OOOO29335 

-0. OOII06 3 8 

.  OO0014881A 

0.273781; 

20 

-0.0039U83 

0.00002361a 

-0.000815U9 

.  OOOOlt331a 

0.136161 

Total 

-1.3913080 

-0. 0010717 ho 

+    .05117920 

0.0000732615 

06 .02^260 

TABLE  VI 
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PRODUCTS  AND  SUMS  OF  BASIC  VAPtlAB] 


FOP  DS3     -   CKATIOH  OP  CT 


*L 


Han  5301 

Point 
Number 

-r<7  (dy) 

j£-  y(ciy) 

°t       dy 

2 

(dy)2 

ll 

0.1X061U 

0.01H2916 

0.002U6957 

3.52633 

.     $ 

0.155313 

0.0350U60 

0.00211369 

7.0U9O3 

6 

0.201780 

O.O5O3380 

0.0007U2736 

7.0U903 

7 

0,221693 

0.065814:0 

-0.000529673 

7.01*903 

0 

0.233375 

O.O012U3O 

-?.00131.1423 

7. 01:903 

9 

0.2236H; 

0.0888060 

-0.00197292 

6.05160 

10 

0.212933 

0.0921816 

-0  .002332141* 

l*.93l02 

.  11 

0.198909 

0.096U313 

-0.00260625 

U.3U723 

12 

0.176299 

0.093873H 

-0.002611492 

3.IJ-510 

13 

0.156U56 

0.039514$ 

-0.002571*80 

2.68960 

lit 

0.129267 

0.0735900 

-  .j;)23^75U 

1.33603 

15 

0.100262 

.     0.07217U6 

-:. 00226538 

l.iil3U8 

16 

0.06101:70 

-O.OOI9616I 

0.938961 

17 

0.069^331 

G. 0529893 

-O.OOI7506I 

0.670761 

18 

0.0569679 

o.ol44#?3 

-O.201i;.7620 

0.1i.502lil 

19 

0 •0^85135 

: .   39217^ 

-0.00131306 

0.337561 

20 

.  335720 

0.O3O77UO 

-0  .001008141 

0.193916 

Total 

2.1:223750 

1.086381i0 

-0.02075300 

59.059300 

BLE  VII  E 


:>9 


PRODUCE  I  SIMS  OF  BASIC  VARIABLES 

FOR  "TAIL  OF?"  OOUSIDE  LA.TIC 


Point 

Huriber 


Run  5301 


y/yd.  t/r  vry  yd  t/r 


(/yd  t/r  )• 

C/CWd  t/r) 


(/yd  t/r  ): 


h 

0.000^60077 

O.OO805O77 

0.00000139216 

0.0000236780 

5 

0,00201119 

0.026361-1* 

0. 00000357 5Uh 

0.000101321* 

6 

0.00520651 

0.01*73322 

0. 00000701*739 

0.000323552 

7 

0.0109806 

0.0778039 

O.OOOOlH-692 

0.000352990 

0 

0.0200673 

0. 115911* 

0.000015931*9 

0.001371*51* 

9 

0.0330517 

0.150723 

O.OOOOI96617 

0.00361-61*7 

10 

0.0$010li9 

0.182653 

0.0000211135 

0.0061*2333 

11 

0.0717637 

0.219186 

0.0000191117 

0.0101*738 

12 

0.0977098 

0.2H2716 

0.0000123288 

0.0161211* 

13 

0.127003 

0.262953 

-0.000000733972 

0.0236298 

12. 

0.159931* 

0.261085 

-0.0000213926 

0.033201*1 

15 

0.191*775 

0.271112 

-0.0000510325 

0. ol*5ioio 

16 

0.23251* 

0.258608 

-0.0000907200 

0.0591*092 

17 

0.27O605 

0.252760 

-O.OOOlla297 

0.0763035 

18 

0.311095 

0.233361 

-0. 00020331a 

0.0958769 

19 

0.352711 

0.23U671 

-0.000277101 

0 •118226 

20 

0.395932 

0.2073^3 

-0.000368757 

0.11*3363 

Total 

2.336101 

3.058936 

-0  .00101*3 

0.63U975 

LE    VIII  E 


.  .OJUCT,,  AND  wlLl:  OF  BALIC  VARIABLES 
FOR  "TAIL  OFF"  GO:.  ATI01I3 


Run  b'301 


60 


Point 

Number 


yw 


?f%  d  t/r 


y2 


1» 

0.209355 

0.0300329301 

0.0132480 

5 

0.523306 

O.OOO0709690 

0.039920 

6 

0.753782 

0.000111679 

0*0825068 

7 

1.001573 

0.00011>761t3U 

0*110353 

8 

1*214088 

0.000171121 

0.21I4823 

9 

1.36616 

O.0OO17321U 

0.299531 

10 

1.U2U76 

0.0O0l6L<732 

0.290833 

11 

1.50109 

0.000130336 

O.U91U71 

12 

1.U7107 

O.00OO7U7233 

0*592207 

13 

1.1:1330 

-O.OOOOOI421362 

0.632606  ' 

Hi 

1*25757 

-0. 00010  30I1.2 

0.770357 

15 

1.17083 

-0.000220391 

0*810161+ 

16 

1.01225 

-3.000355097 

0.910211 

17 

0.89633!; 

-0.000501065 

0.959616 

13 

0.773U13 

-0.000659736 

1.3G9U2 

19 

0.700109 

-O.OOO826692 

1.05227 

20 

0.57263 

-0.0010135 

1.09370 

Total 

17.293it67 

-0.002601 

9.535292 

61 


TABLE 


jjl  h 


PHODDGTS  AIID  SUMS  0?  BAolC  VA  SABLES 
FOR  !'TAIL  OF?1'  CONSIDERATIONS 

Ran  >301 


Point 

Number 


w 


Aj^  d  t/r     (  /cnt  dt/r  ): 


u 

0*000520307 

0.0000000313532 

.5 

0.000930322 

0.000000126167 

6 

0,00102707 

0.000000151165 

7 

O.OOlOiiolii. 

0.00000015U213 

0 

0. O00933l4l;l 

0.000000136309 

9 

0.000312698 

0.000000106015 

10 

0.000600517 

0.000000069U32 

11 

0.000399762 

0.0000000 3U8569 

12 

0.000135616 

0.000O0OO09l|281; 

13 

-0 •0000087 2ho6 

0.00000000002601 

lit 

-0.000168211 

0.000000013783 

15 

-0.000306767 

0. 000000057 Ikh 

16 

-0.00039U90U 

0.000000138533 

17 

-0.000U67702 

0.000000261271; 

13 

-0.000505523 

O.OOOOOOU31255 

19 

-0.000550027 

0.0OO00O6it9U75 

20 

-O.000 53 3303 

0.0000009l431;3l 

Total 

0.0035766 

0. 00000 339U1 

TABLE  ::  E 


L>< 


Point 
Number 


iDUCTO  A! ID  SUMS  OF  BA^IC  VAHIAB] 

Fort  j:;t:  ^to:;  0/  lumped  constants 

Run  5301 
nz  ^/nzd  t/r  (  /r^d  t/r  )</£d  t/r  )  jj^d     t/ 


2 

0.00001003901; 

-0.3000002^3116 

-0.00000770381* 

!| 

0.000136517 

V-0.000001509UO 

+  0.000011183692 

6 

O.OOU07130 

0.000025308^ 

0.00011*91*91 

8 

0.0187U96 

0.0000979277 

0.000U22363 

.10 

0.0U9U262 

0.  0O02hl«OU7 

0.000823513 

12 

O.0?36ol; 

0.3-001*86670 

0.00133615 

1U 

0.16U261; 

0.000839U91- 

0.00200l*2l* 

16 

0.21,11.31 

0.00131379 

0.00272315 

13 

0.32l|695 

0.00190013 

0.0031*0961 

20 

0.1al;297 

0.002623Q8 

0.001*13133 

Total 

1.3l60li2 

0.0075^050 

0.0150570 

TABLE  XI  E 


63 


PHODOCTS  AHD  SUMS  OF  BA^IC  VA3XABLES 
FOR  DETECT,  IATIOIT  OF  LUMPED  COHETA1TTS 

Run  5301 


Point 
Number' 


hdr^Jn^d   t/r       (  r^d  t/<r  )2        r^  /T  d  t/j 


2 

^0.0000100036 

0.00000050694U 

-0.  0000050 37  hk 

b 

*O.OOOOli5l2i>3 

0.00O00203913 

-/-0. 00010101*9 

6 

0.0005U5305 

0. 00020701^5 

O.OOOi-97797 

0 

0.001i;25l3 

0.001581*63 

0.00115866 

10 

0.00232095 

0.0059U56U 

0.00202377 

12 

0.0032lj655 

0.0156520 

0.00307699 

Ik 

0.003l|l|l55 

0.03300178 

O.OOla76U8 

16 

0.0032697U 

0.0600779 

0.0052796 

10 

0.00307177 

O.O970i;38 

0.00633130 

20 

0.00137963 

0.1U7155 

0.00738720 

Total 

0.0192358 

0.361U91 

0.03003331 

TABLE  XII 


I     )DUCTS  AIID  SUMS  OF  BASIC  VA:2ABLEo 
5    I  DETES12HATIOH  07  LUMPED  GOIJ^TAilTG 


Ihm  5301 


6U 


Point 

Number 


flu  :<n 


n^hdrij 


"a 


2 

-O.OOC 

?3315 

-0.000199089 

O.OOO2OO709 

I* 

■/•o.ooc 

>323U 

+  0.00302096 

0.00913936 

6 

o.ooi 

t037 

0.0107257 

0.0800390 

8 

o.ool 

?731 

0.0168618 

0.2218U1 

10 

o.oo£ 

^88 

0. 019291a 

O.iaoOGl 

12 

O.OOC 

1.730 

Q.Q2QZ2&& 

0.625631 

u* 

0.005 

fll2 

0.0171218 

0.317216 

16 

0.01C 

Oh 

0.0131399 

0.970225 

18 

0.013 

bi2 

0.0101932 

1.077U1*) 

20 

o.oii 

'20 

0.00^2920 

I.l662j0 

Total 

0.068 

'21 

0.115977 

5.379113 

table  nn  s 


65 


FOR  DETEiC  TIJATIOH  0?  LUMPED  COiToTAIJTo 


Point 

Number 


Hun  ^301 


Sfid  t/r 


hdn, 


s /dd  t/r 


c  y^ci  t/v)2 


2 

0.00000  38)4651 

O.OOOOOU99U73 

0.000000126300 

u 

0,0000109022 

0.000033U012 

0.00000111725 

6 

0.0000102760 

0.0000666661 

0. 00000  309J408 

8 

0.0000261006 

0.0000600680 

0.00000605160 

10 

O.OOOO330O22 

0.0000952665 

0.00000017225 

12 

O.OOOOUl5k52 

0.0001009^55 

0.0000151321 

lit 

0.0000509586 

0.0000875020 

0.0000213liljl| 

16 

0 .oooo595k96 

0.0OO071502U 

0.0000207296 

10 

0.000066U90 

0.0000599020 

0.000037210 

20 

o.oooo7U556o 

0.0000335160 

O.oooolj.67056 

Total 

0.000386107 

0.0006L11765 

0.000169603 

TABLE  XIV  E 


PxiODUCTo  AIID  SUMS  0?  BASIC  VARIABLE.; 
POR  DBXSRUINATZON  0?  LU1IPED  CONSTANTS 


Run  5301 


Point 

iiumber 


S  (lidxig) 


2 

0.000152021 

0.000117072 

u 

0.00032832U 

0.000107952 

6 

0.000393781 

0.000107952 

8 

0.000379833 

0.000112572 

10 

0.000321^60 

0.OO011U062 

12 

0.0OO277H6 

O.OOOIH1O62 

lit 

0.000203908 

0.000121661 

16 

0.00011$  207 

0.000123U321 

13 

0.000107023 

0.000113310 

20 

0.000053)410 

0.000113810 

Total 

0. 00237011a 

0.001156385 
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TABLE  XV 

BAoIC  VARIABLES 
Hun  5068 


Point       Time 
ITunber     (sec.) 


**(z)      *s     ^rad*) 


.# 


df 


2 
3 
2i 

6 

7 

8 

9 
10 
11 
12 
13 
U* 
15 
16 

17 

10 

19 
20 


0.10 
0.1> 
0.20 
0.2^ 
0.30 
0.35 

o.Uo 

o.U5 
o.5o 

0.55 

0.60 

o.65 

0.70 

0.75 
o.8o 

0.35 

0.90 

0.95 
1.00 


-O.Ollili; 

-f  0.0319 
.  0.0921 

0.173U 

0.21-63 

0.322 

0.396 

0.U78 

0.1-97 

0.U9U 

0.U61+ 

0.U2G 

0.390 

0.350 

0.298 

0.21-8 

0.209 

0.159 

0.122 


1.3U0 

1.132 

1.953 

1.953 

2.160 

2.160 

2.160 

0.519 

0. 

-O.308 

-0.772 

-1.231; 

-1.23U 

-1.23U 

-1.23U 

-1.231; 

-1.338 

-1.338 

-0.722 


-0.0262 
-0.033U 
-O.0U9U 
-0.0625 
-O.0639 
-0.0639 
-0.0639 
-0.0523 
-0.0392 
-0.0233 
-0.0138 
-0.00372 
O.OO36U 
0.00651; 
0.00726 
0.00372 
O.OO872 
0.00372 
0.00372 


-0.0271 
-O.031U 
-0.01-35 
-0.05H; 
-O.OU32 

-0.01-35 

-0.0336 

-0.0220 

-0.0071; 

O.OO83 

0.0159 

0.0106 

0.0235 

0.0289 

0.0263 

O.02I46 

0.0221 

0.0189 

0.0165 


-0.02U35 
-0.02020 
0 

0.01728 
-0.00U33 

0 

-0.26200 

7^0.0201-5 

+0.02000 

0 

0.001-31 

0 
0 
0 
0 
0 
0 
0 
0 


nz 


-rcj 


Oo 


le  m  s 

BAOIC  VARIABLES 

Run  5063 

Lnt 

Ifm  iber 

Time 

(sOC.) 

r 

CL 

2           ^ 

1 

dy 

2 

0.10 

-0.00397 

-0. 00025 

-0.001158 

-1.I48 

3 

0.1$ 

-0.01*357 

-O.00278 

-O.D007U6 

-1.016 

h 

0.20 

-0.009938 

-0.00634 

v-o.  000U65 

-2.090 

5 

0.25 

-0.17316 

-0 .01103 

O.OO0O15U5 

-2.090 

6 

0.30 

-o.°.l>5l<7 

-0.01568 

0.  )00o5305 

-2.090 

7 

0.35 

-0.31579 

-0.02015 

0.000398 

-2.090 

3 

o.Uo 

-0.37731 

-0.02U10 

0.001193 

-2.090 

9 

o.h$ 

-O.U5205 

-0.02390 

0.001660 

-1.600 

10 

o.5o 

-0.47100 

-O.O30IO 

0.1662 

-0.262 

11 

0.55  . 

-0.47130' 

-0.03015 

0.001I1.20 

y-0.262 

12 

o.6o 

-0.  W  '1-20 

-0.02335 

0.001265 

0.386 

13 

o.6$ 

-O.iOlOii 

-0.02630 

0.001032 

1.091 

in 

0.70 

-0.3756U 

-O.024OO 

0.000336 

1.211 

15 

0.75 

-O.3388O 

-0.02165 

0.000716 

1.211 

16 

O.oO 

-0.28902 

-  1  .OI8J46 

0.0005U2 

1.211 

17 

o.85 

-0.21j236 

-O.01&9 

0.000361 

1.211 

18 

0.90 

-0.20550 

-0.01311 

0.000224 

1.211 

19 

o.?5 

-0.15679 

-0.01002 

0.000161 

1.211 

20 

1.00 

-0.12016 

-J. 00770 

0.000104 

1.211 

'H 


-  nr 


TABLS  XVII  E 

BASIC  VA  SIABI3S 

iiun  >068 


69 


Point      Tune 
IJiiraber     (sec.) 


Jyd  ify 


w 


fcH  d  t/r  /? 


dt 


2 

0.10 

-0.0002065 

-1.092U 

-0.00007672 

-0.00126 

3 

0.15 

-0.0009735 

-1.7302 

-0  .00011*312 

-0  .00221* 

U 

0.20 

-0.0033' 35 

-2.112a 

-O.0O013UO2 

-0.00372 

5 

0.25 

-0.0031235 

-2.0907 

-0.00019932 

-0.00565 

6 

0.30 

-0.0152535 

-2.0928 

-0.00019707 

-0.00793 

7 

0.35 

-0.02U3035 

-2.1091 

-0.00013331 

-0.01023 

3 

o.Uo 

-0.0367235 

-2.1U73 

-0.00013171 

-0.01236 

9 

0.16 

-0.0511235 

-1.6799 

-0.00003191 

-0.01U28 

10 

o.5o 

-0.0670235 

-0.310.5 

+0 .00008609 

-0.01590 

11 

0o> 

-0.0332235 

V-0.1933 

0.00019209 

-0.01696 

12 

0.60 

-0.0991235 

0.3253 

0.00028219 

-0.01752 

13 

o.65 

-0.1138735 

l.oiaU 

0.00036029 

-0.01733 

H; 

0.70 

-0.1273735 

1.1681; 

0.0001-2629 

-0.01799 

15 

0.75 

-0.139ll035 

1.1766 

0.000lt8069 

-0.01781 

16 

0.80 

-0.1500035 

1.1850 

0.00052269 

-0.01758 

17 

o.05 

-0.1591735 

1.1937 

0.00055399 

-0.01725 

18 

0.90 

-0.1669235 

1.2002 

0.00057H6U 

-0.01697 

19 

0.95 

-c. 17 31335 

1.2033 

0.00053339 

-0.01665 

20 

1.00 

-0.1700U35 

1.2060 

0.0005973U 

-c .01639 

*    -     djT   -    (CL01   /    20J,    )     n 

Xi 


tabus  r/in  e 


70 


BASIC  VAillABLBS 

Kun  5068 


Point      Time  /(rad.) 

1  Junker     (sec.) 


/Jd  t/r       hdn2 


ykgd  t/-, 


2 

0.10 

0.00926 

0.000314U 

0.01915 

-  .000757 

u 

0.20 

0.00790 

0.000931 

0.027U 

■f  0.001303 

6 

0.30 

0.00837 

0.0011*92 

0.0308 

0.0133U3 

0 

O.iiO 

0.00185 

0.001972 

0.0283 

0.0355U3 

10 

0.50 

0.00023 

0.001957 • 

0.0021 

0.067793 

11 

0.55 

0.00028 

0.001966 

-o.ooia 

0.031*693 

12 

0.60 

0.00023 

0.001966 

-0.0113 

0. 1011*1*3 

13 

0.65 

0 

0.001966 

-0.0171 

0.1167U3 

li; 

0.70 

0 

0.001966 

-0.0171 

0.1307l;3 

15 

0.75 

0 

0.001966 

-0.0171 

O.li.3253 

16 

0.80 

0 

0.001966 

-0.0253 

0.153993 

17 

0.35 

0 

0.001966 

-0.0171 

0.16&53 

18 

o.?o 

0 

0.001966 

-0.0173 

0.171393 

19 

0.95 

0 

0.001966 

-0.0151 

0.177863 

20 

1.00 

0 

0.001966 

-0.0113 

0.182603 

TABLE  XIX  E 


71 


PRODUCTS  AND  SlttE  OF  BASIC  VARIABLE 
FOR  DET  ^TilATION  OF  CT  J 


Run  5068 


Point 
Number 


vdr^ 


vd^ 


(dXJCdna) 


(d<T)' 


(dnjs)' 


2 
3 
!» 

6 

7 

8 

9 

11 

12 

13 

12* 

15 

16 

17 

18 

19 

20 


-0.036311; 

-0.03&'UU8 

-O.C8U9555 

-0,10038U2 

-0.10iO.120 

-0.0939600 

-0.0333760 

-0.C11U200 

-O.G02556U 

-0 .0122718 

-O.C22952U 

-0.0351690 

-0.0356626 

-0.032l£U2 

-O.03O356U 

-O.O295698 

-0.0252832 

-0.0119130 


O.0O0673U35 

0.00063it230 

0 
-0.000888192 
•fO.OOO2O07O6 

0 

0.01011320 

-o.oooii5iooo 

0 

+0.000063529 
0 
0 
0 
0 
0 
0 
0 
0 


-0.033299 
-0.022866U 

0 
f0.0337U73U 
-0.00935280 

0 

-0.5659200 
+ 0.01061355 

0 

-0,00332732 


0.0006160 

O.OOOliOSO 

0 

0.0002982 

0.0000188 

0 

0.0636000 

O.OOOU185 

0 

o.ooooi355 

0 
0 
0 
0 
0 
0 
0 
0 


1.795600 
1.23LU2U 
3.31U209 
3.81U209 
U.665600 
U.665600 
U.665600 
0.269361 
0.09U86U 
0.59593U 
1.522756 
1.522756 
1.522756 
1.522756 
1.522756 
1.790210; 
1.79021.J; 
0.52128U 


Total        -O.7882638      0.010359133      -0.590hOJil3        0.07037805  37-373003 


£  XX 


"product:;    a;: a  suib  of  ba^ic  va^iable^ 

FOR  DEXSRHENATION  OF  CL, 


Run  5068 


j0i 


72 


Point 
Number 

T-^(dy) 

2       * 

0L    y(dy) 

Mr)* 

2 

0.0300776 

0.00132933 

0.000290Ui; 

1.31790 

3 

o.o6o6£m- 

O.0O135U7U 

0.0050U811.3 

3.29736 

It 

0.1032I16 

-0.000971350 

0.0132506 

U. 36310 

5 

0.130625 

-0.0000322905 

0.0231572 

U.36610 

6 

0.133551 

-0.00011007U 

0.0327712 

h* 36310 

7 

0.133551 

-0.000331820 

0.0U21135 

ii.  36310 

8 

0.133551 

-0.002^9337 

0.0503690 

U.36310 

9 

0.0836800 

-0.00265600 

0.0U62i±0Q    • 

2.5600 

11 

-0 .00667988 

+ 0.000371330 

-0.00733U23 

0.063382 

12 

-0.0122268 

0.00112079 

-0.0251181 

0.73U996 

13 

-0.00951352 

0.00112591 

-0.0236933 

1.190230 

m 

*O.OOUU0199 

0.00107295 

-O.O29O61;0 

1.1-6652 

15 

•  0.0079199k 

0.000367076 

-0.0262132 

1.1:6652 

16  , 

0.00379136 

0.000656362 

-0.0223551 

1.1*6652 

17 

0.0105599 

0.0OOU37171 

-0.018758U 

I.I16652 

18 

0.0105599 

0.000293062 

-0.0158762 

1.U6652 

19 

0.0105599 

0.00019U729 

-0.01213U2 

1.U6652 

20 

0.0105599 

0.000126307 

-0.00932U70 

1.1-6652 

Total 

•  O.OW469 

0.00135360 

o.oiyOiitO 

la.325553 
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tablb  :sa  3 


PRODUCTS  AND  'M£j  OS  BASIC  VARIABLES 
SOU  "TAIL  OB71'  CONSIDERATIONS 


Point 

Iluriber 


Rim  5068 


^  t/r 


-y^d  t/7 


(jjd  t/-r){fcntd  t/r) 


2 

O.OOOOOO8190O5 

0.000225581 

0.00000001531*27 

3 

o.ooooi<2ia5U 

0.00173302 

0.000000139327 

Ii 

0.000332260 

0.0070626U 

0.000000615271 

5 

0.0011.066^ 

0. 016901a 

0.00000162321; 

0 

0.00371^26 

0.031922U 

0.00000300601 

7 

0.00783270 

0.0523131 

O.0OO0OU54673 

8 

0 .0138261 

0.0788561* 

O.OOOOOU83635 

9 

0.02311014 

0.0058821. 

0.00000163135 

11 

0.039261;8 

-O.OI0OO7I 

-0  .0000159861; 

12 

O.OUU0307 

-O.OO1OOI7 

-0.0000279717 

13 

0.0)468977 

-O.110508 

-0.00001*10275 

Ui 

O.OU78I*66 

-0.11*8823 

-O.OOOO5U290O 

15 

0.01*72299 

-0.161*022 

-0.0000670099 

16 

O.OU335140 

-0.177751* 

-0. 00007  31*o53 

17 

0.0385773 

-0. 180002 

-0.0000881305 

10 

0.03U3019 

-0.2003i{6 

-0.0000959209 

19 

0.0271526 

-0.203388 

-0.000101899 

20 

0.0213937 

-0.211*720 

-0.000106353 

Total 

0.1ii0375 

-1.235552 

-0.000660637 

71* 


TABLE  .:::;■; 


PRODUCTS  AMD  STBS  OF  BAoIC  VAilIABLSO 
FOR  "TAIL  OFF"  COHSIDS3UTIONS 


Run  p06G 


Point 
Number 


(fjd  t/r  ): 


yvr 


y/cnt  d  tA 


2 

0.00000001*261*23 

0.002*33683 

0.000000301*578 

3 

0. 00000091*7702 

0.0775633 

0. 000006  23571* 

u 

0.0000111790 

0.209915 

0.0000132870 

5 

0.0000659913 

0.362029 

0.000031*6001* 

6 

0.000232669 

0.513715 

O.OOOOi*837ii6 

7 

O.OG06l52ll> 

0.666033 

0.0O0O578375 

3 

0 ,00131*862 

0.810198 

0.00001*96955 

9 

O.OO26136I 

0.759399 

0.000011*1*21*9 

11 

0.00692613' 

-0.0911939 

-0.0000906231 

12 

0.0093251*7 

-0.366576 

-0.00012531*9 

13 

0.0129672 

-0.1*23390 

-0.00011*8382 

ll* 

0. 01622L.0 

-0.1-33898 

-0.000160132 

13' 

0.0194333 

-0.393632 

-0.000162853 

16 

0.0225105 

-0.31*21*89 

-3.000151068 

17 

0.0253362 

-0.239300 

-0. 000131*265 

18 

0.0278635 

-0.2U66I*1 

-0 .000113036 

19 

0.0299925 

-0.183656 

-0.0000922507 

20 

0.0316995 

-0.11*1*913 

-0.0000717761* 

Total 

0.207657 

-^0.1^6995 

-0.001021*98 

:'. -,'  u  i    AAXJ.X 


products  a:;d  sums  of  ba^ic  variabi 
for  "tail  0j7!<  coiish 


Point    ' 
Ilimiber 


y2 


I&n  5068 


w/Cj^  d  t/r  (/cnt  d  t/r  ): 


2 

0 •00001^7609 

0.0000838089 

0.00000000588596 

3 

0.0013933U 

0.0002^732 

0.000000020)1833 

u 

0.00937^39 

0.000333715 

0. 00000003336 3h 

5 

0.0 2 990 37 

O.OOOl<17772 

0.0000000399230 

6 

0.0602550 

O.OOOla2li26 

O.OOOOOOO3GG366 

7 

0.0997233 

0.000386619 

O.OOOOOOO336026 

3 

O.U  I2363 

0.000232821 

O.0OO0000173U75 

9 

0.201-3U9 

0.0000536056 

0.00000000101825 

11 

0,222595 

0.0000371310 

0.0000000363936 

12 

0.1973H; 

0.300232377 

0.0000000796312 

13 

0.169612 

c. -00375206 

0.000000129309 

Hi 

O.U4II05 

0.000Ji98077 

O.OOOOOOI81723 

15 

o.ni.735 

0.000565580 

0.000000231063 

16 

0.0835326 

0.000619388 

0.000000273205 

17 

0.058733U 

0.000661287 

O.OOOOOO306905 

13 

0  .01*22282 

0.000689700 

0.000000330211 

19 

0.02U5815 

0.00070799O 

0.0000003U6203 

20 

0.01iili38U 

0.000720392 

0.000000356315 

Total 

1.617391; 

0.007339376 

0.000002^631:3 

TABLE  XHV  E 
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Point 
Lkinbcr 


PRODUCTS     AND  STUB  OP  BASIC  VAPIADLEo 

POP  DPPP.il IIIATIOII  OP  L01SPSD  CGNSTAHTS 

Run  5063 
(d£    and  £  assumed  to  »  0) 


nz/n^l  t/r  ( /nzd  t/r>( // d  t/r)       hdx^J^d  t/T 


2 

0.00001070U 

-0.00O00O260l|03 

-0.000013i;397 

1* 

0.000120006 

+ 0.00000121309 

-/•  O.OOOO33378O 

6 

0.00323638 

0.0000199078 

0.000333611 

8 

0.01it07>0 

0.0000700903 

0.00102136 

10 

0.0336931 

0.000132671 

0 

12 

0.0U70696 

0.000199U37 

-O.OOlOiiOol 

Hi 

0.05092Ui 

O.0002570hl 

-0.0021J45U9 

16 

0.0U  53129 

0.0OO3O275O 

-0.00252703 

13 

0.03582U 

0.000336959 

-0.00305030 

20 

0.0222IP.0 

0.000353997 

-0.00175299 

Total 

0.25305U 

0.00167391 

-0.00909126 

TABL2  XXV  E 
PRODUCTS  AND  SOI &IC  VARIABLES 

for  02::     :     Tim  of  lumped  co:;v;ta. 

(<lS    and   <T  assumed  to  ■  0) 


Point 
iber 


(JV  t/ 


■T  / 


% 


/7  d  0/7 


v*3^ 


2 

0.000000573049 

-0.00000486416 

-0.000251975 

u 

0.000001697^1 

^0.00003571^1 

+  '.00239U60 

6 

O.OOOI73036 

0.000367^0 

0.00708113 

8 

O.OOI2633O 

0.000780912 

0.011385 

10 

o9ooh$9&9 

0.000972S29 

0 

12 

0.0102907 

0.000912224 

-O.OO476064 

Hi 

0.0170937 

0.000765757 

-O.OO639170 

16 

0.0237138 

O.Ojo531<G85 

-0.00408198 

18 

0.0293756 

'  O.OOOI1IO094 

-0.00372020 

20 

0.0333439 

0.000239U59 

-0.00116928 

Total 

0.119356 

0.00511512 

-0.0003150U5 

(o 


LE  XXVI  E 

PiiODUCTo     AID  ...UTL  OF  BASIC 

Foh  det  :oh  07  lumped  go: 

Pain  5068 
(di"    and  <f  assumed  to  =  u) 


Point 
I  lumber 


*V 


hdng  yVd  t/r  (  /i"  d  t/r  ): 


2 

0.0001999l|0 

0.00000613008 

O.00OO0011G336 

U 

0.008U821O. 

O.00002i;2060 

0.000000366761 

6 

0.0606637 

O.OOOOU23950 

0.00000222606 

3 

0.156816 

O.OOOC566950 

O.OOOOO333G78 

10 

0.2U7009 

0 

0.00000332935 

12 

0,215296 

-0.0000201712 

0.00000386516 

lit 

0.151710 

-0.0000322621 

0.00000386516 

16 

0.0335063 

-0.0000322621 

0.00000336516 

13 

O.OU368I 

-0.00003U99U8 

0.00000386516 

20 

0.01U8352 

-0.0000138736 

0.00000386516 

Total 

0.987200 

-O.0OOO138736 

0.00000386516 

JU  21  62 
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Thesis 
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Thesis  18194 

Vt69  les 

tA  &  -_  of   the  equations 

of  motion  method  for  the  evalua- 
tion of  the  longitudinal  dynamic 

bilit  derivatives  of  the  F-8CA 
fror,i  transient  response  data. 
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